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^2 . Abstract. Classical chaos refers to the property of trajectories to diverge 

^ I exponentially as time t — > oo. It is characterized by a positive Lyapunov 

Q ' exponent. 

• There are many different descriptions of quantum chaos. The one re- 

, i lated to the notion of generalized (quantum) Lyapunov exponent is based 

^ I either on qualitative physical considerations or on the so-called symplectic 

^ ' tomography map ^ |2 . 

. The purpose of this note is to show how the definition of quantum Lya- 

l/^ I punov exponent naturally arises in the framework of the Moyal phase space 

■ formulation of quantum mechanics [Sj , and is based on the notions of quan- 

^ . tum trajectories and the family of quantizers [1]. The role of the Heisenberg 

I uncertainty principle in the statement of the criteria for quantum chaos is 

p ' made explicit. 

C ' 1. Introduction. 

■ 

^ ; 

qh, Irregular behaviour of classical dynamical systems arising from deter- 

^ ' ministic time evolution without any external randomness and stochasticity 

- the so called deterministic chaos - manifests itself as an extremely sensitive 
dependence on the initial conditions, which makes unstable the long-time 
I prediction of the dynamics. 

In such a system, a positive Lyapunov exponent is a quantitative measure 
of the infinite time exponential separation of neighbouring orbits. 

In detail, let the system have a form x = F{x), where x = {q,p) is 
(i-dimensional vector from the system's phase space. Denote by x{t,xo) 
its solution with initial point at xq and t G (0, oo). Then the Lyapunov 
exponent is given by 

, 1, llfaWII 

Here || - H represents the d-dimensional Euclidian norm, and ||(5x(0)|| is initial 
infinitesimal deviation from xq, \\6x{t)\\ = \\x{t,xo -\- 5x{0)) — x(t,XQ)\\ is 
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deviation from x{t,XQ) at time t. In the limit ||5x(0)|| — > one obtains 

A^, = lim ]■ In [{v ■ Vxo)x{t, xq)] , (1) 

where v G R'^ is a unit vector in the direction of the initial displacement, 
5x{0), and Vxq is gradient w.r.t. the initial point xq. 

To extend this notion to quantum mechanics it is natural to use its 
phase space formulation where quantum observables on the Hilbert space 
Ti. = L^(M"',C) are represented by functions on the phase space, called 
symbols. Consider Weyl symbol A{x), obtained from the operator A s.t. 
M{q') = / {q'\A\q"){q"\'il^)dq", V G by formula 

A{x) = 2" y e^'P'/^iq - s\A\q + s)ds = [A]^{x) . (2) 

The dimension of the phase space is even, d = 2n. 
In this case quantum mechanical mean value 

= Tr Ap= I {q'\A\q"){q"\p\q') dq'dq" (3) 

can be written in a form analogous to classical statistical mechanics 

{A)p = j A{x) W{x) dx , (4) 

which gives mean value of A in the state p{x) = [p]vi{x). Here W is the 
Wigner distribution function W{x) = h~^p{x). 

One drawback of this approach, relative to the classical statistical den- 
sity, is that the W{x) is not everywhere nonnegative so is not a conventional 
probability density. 

To circumvent this difficulty another representation can be considered 
in which a new density is defined as the Radon transform 5 of the Wigner 
function W{x). It takes only non-negative values and becomes straightfor- 
ward analog of the classical statistical density. In it was shown that the 
tomographic representation of quantum mechanics based on Radon trans- 
form is an alternative to the Weyl-Wigner formalism. 

Due to their similar nature in quantum and classical cases the tomo- 
graphic distributions were used in to define a quantum Lyapunov expo- 
nent. In detail, let U{t) = exp [—itH/h] be the Schrodinger unitary evolution 
operator, and X{t) = U^{t) x U{t) be time evolution of the quantum coor- 
dinate operators x = {p,q)- The formula for quantum Lyapunov exponent 
turned out to be 

A, = ^lim iln||(X(t))ps||, (5) 
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where average is taken with respect to a special initial singular density p^, 
the kernel of which has the form 

WW) = ((»,v) + 'J;li^pn.'J , (,„ _ ^si-j . 

Here parameters {qo,Po) correspond to the initial point xq in the phase space, 
and vector v = {vi,V2) ^ defines the direction of the initial deviation from 
Xq. The authors of [2] stress special role of the tomographic distributions 
for obtaining these formulas, in particular the one for the initial density 
In similar formula is obtained from qualitative considerations but the 
choice of the initial density appears rather ad hoc in the quantum mechanical 
setting. As it is shown in ^[2] , quantum Lyapunov exponent helps to 
classify different types of quantum complexity. There are examples where 
exponential rate of growth for the trace {X{t))pv = TrX {t)p^ of position and 
momentum observables starting from the singular initial density matrix © 
was found in quantum mechanics. In many cases when quantum mechanics 
has damping effect on the classical chaos and the rate of growth for the trace 
is milder than exponential, the notion of quantum sensitive dependence was 
used instead. 

In this paper we derive representations of the quantum Lyapunov index 
from the Weyl-Stratonovich quantizer ^ Using notion of quantum tra- 
jectory 0, i.e. the symbol of operator X{t) , we rewrite the formula for 
quantum Lyapunov exponent in the form identical to the classical definition 

(see formula (HU), replacing classical trajectory x{t,XQ) with quantum 
trajectory X{xo,t;h). In this form it becomes transparent that the defini- 
tion respects the correspondence principle: in the limit ^ — > definition of 
quantum Lyapunov exponent transforms into the classical one. In contrast 
to the classical trajectory, the Heisenberg uncertainty principle prevents the 
quantum trajectory to be interpreted as a measurable physical value despite 
the fact that it can be viewed as a limit of a sequence of quantum means. 

2. Quantum Phase Space. Weyl Stratonovich quantizer. 

Development of phase space formulation of quantum mechanics has a 
long history, but it still is of interest due to extensive study of possible 
generalizations to the case of non-Abelian gauge theory with support on a 
Riemannian manifold. 

In 1932 Wigner introduced his quasi-probability distribution associated 
with the wave function ^{q) 

W{x) = (l) [ e'^'P'^^iP+iq- s)ij{q + s)ds, h = 2TTh. (7) 
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Further developing Groenewold's ideas published in 1946, Moyal gave 
in his paper of 1949 statistical interpretation of the Wigner's formula as a 
Fourier inverse of the expectation value of the Heisenberg translation oper- 
ator 

T{y) = exp{—2iJy ■ x/h) , T{x)^ xT{x) = x — 2x1 . (8) 

namely 



where J denotes the Poisson matrix 



He also showed that the 



/ 
-/ 

Wigner rule ((TJ of getting the phase space function from the operator 
(and A) is inverse to the Weyl quantization rule, 

{q'\A\q") = h-- [ e^vi^-'i'^I^A {^^.^ dp. (9) 



A significant step further was the introduction |Sj of the family 
of unitary operators A{x) labelled by points of the phase space. These 
operators define both quantization © and dequantization ((2} rules, which 
establishes a unitary isomorphism between symbols and operators. 



A{x) = 2" Tr A(x) A, A = j A{x) A(x) d*x . 



(10) 



For this reason the operators /S.{x) are called quantizers. Above d*x = 
{'Kh)~'^dx denotes a dimensionless phase space measure. 

The fundamental nature of quantizers also reveals itself in the fact that 
they define the noncommutative product * for the phase space functions 

{A*B){x) = j A{y)B{z)K{x,y,z)dydz, (11) 



2*^ 

K{x,y,z) =^-^T¥A(x)A(y)A(z). (12) 



Other useful properties W of quantizers are: 

1. A(x) = A(x)+ = A(2;)~^ Thus A(2;)2 = / and ||A(a;)|| = 1 : 

2. TrA(x) = 2-"; 

3. Tr [A(x) A{x')] = {^h/2Y 8{x - x') ; 

4. A(x) = fe'^''' -^y/^T(y)d*y; 
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5. p{x) = 2"(V'|A(x)|V'> ; 

6. A(0)x A(0) = -x; 

7. A(x) = r(x/2)+A(0)r(x/2); 

8. fA{x)d*x = I. 

Many authors (see e.g. ^0]) use quantizers as a fundamental object 
defining the deformation quantization introduced in 1978 by Bayen et al 

nu. 

There is one observation useful for the purpose of our note. In view of 
the 1st equation of (|1()|) Property 3 tells us that the symbol of the quantizer 
is [A(x)]w(x') = (/i/2)'^5(x - x'), and by © that its kernel is 

{q'\A{x)\q") = 2-"e*P{9'-9")/^<5 (^g _ ^-^^ . (13) 

An attractive feature of quantum phase space is computation of the trace 
of an operator and pairs of operators. For an operators A, B with symbols 
A{x),B{x) one has 

Trl = ^ j A{x)dx, (14) 

TtAB = j A* B{x)dx = ^ J A{x)B{x)dx . (15) 

Identity (|14j) follows from (|in|l and Property 2. The removal of the * oper- 
ation in ()15() is a consequence of Properties 3 and 8. 

3. Quantum means and symbols. 

Mean value of a quantum observable given by operator A in a unit nor- 
malized quantum state ■0(g) is 

(V^lllV^) = TVlp, p=m{i^\. (16) 

This is an example of formula © for the pure state density. Its phase space 
form is (j^ which is a special case of the trace identity (|15j) . 

Let us consider a family of Gaussian states localized near q = qq £ i?" 
with width \/^ 

Wh[q; qo,po) = , exp — — + -po{q - go. 
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These states are all unit normalized, ||'0/i('i'O)Po)|| = 1- By formula (|7|) cor- 
responding normalized Wigner function has the form 



Wnix;qo,po) = — — exp 



{q - qo? + {p- po? 



Two important remarks should be made about this function. 

Remark 1. Function Wfi{x] xq) is positive and thus can be interpreted as 
the classical statistical density. 

Remark 2. In view of the following formula 

lim Whix; xq) = 6{x - xq) (17) 

it is evident that the symbol of the quantizer A(x) is proportional to the 
limit of the sequence of Wigner functions Wfi{x; xq) of the localized Gaussian 
states as their width \/^ goes to 0. 

Thus the /i-independent symbol A{xq) of an operator A evaluated at 
point xq of the phase space appears as the mean value calculated with respect 
to the localized Gaussian state with parameters xq = {qo,Po) in the limit 
when the width of the state goes to zero 



lim{'il;n{q;xo)\A\ipfi{q;xo)) = lim / A{x)Wh{x; xq) dx = A{xo). 

Note that this limiting process takes one outside the framework of quan- 
tum expectation values. The norm ||A(a;)|| = 1 implies that for any pure 
state p = |V')(V'I 

2X n 



|p(x)|<2- and \W{x)\< ^-j . (18) 

This means that the quantizer or any singular symbol such as 6{x — xo) does 
not correspond to a pure state. The global bound ((TH|) forces p{x) to be 
distributed in phase space without large peaks and is an evident consequence 
of the uncertainty principle. 

Let quantum operator A have Weyl symbol A{x), and time dependent 
operator A{t) = U+{t)AU{t) = A{X{t)) has Weyl symbol A{x,t;h) such 
that A{x,0;0) = A{x). Consider now the following expression 



{A){t,xo;h,e) = / We{x;xo)A{x,t;h) dx 
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as a function of two small parameters h and e. The following table shows 
meaning of the expression if one or both of the parameters are 0. 



{A){t,xo;h, e) 







e ^ 


classical observable 
A{x{t,xo)) 


symbol A{xQ,t; h) of 
quantum observable A{t) 


e / 


classical statistical mean 
fdyWe{y;xo)A{x{t,y)) 


quantum mechanical 
mean value 



Note that function ioi: e ^ h does not correspond to a pure state 
ip{q), but rather is a density matrix. 

4. Classical and quantum Lyapunov exponents. 

From the table we see that the symbol A{xQ,t;h) is in fact quantum 
analog of the classical value A{x{t,XQ)). This is due to Egorov Theorem 
which states 

lim A{x, t; h) = A{x{t, xq)), (19) 

?i— >0 

where x(t, xq) is classical trajectory. 

Thus the quantum analog of {v-V xa)-^{x{t-, xq)) will be {v-V xo)-^{xQ^t; h). 
From the 1st formula of (|1U() one has 



A{xo,t; h) 



{q'\^{x^)\q"){q"\A{t)\q')dq"dq' ^ 



Here we see that whole information about the point xq of phase space is 
contained in the quantizer kernel {q'\/S.{xQ)\q") . Thus the derivative (w-V^^q) 
will only effect this part of the formula 



{v ■ V,.„)A(xo, t- h) = r j [{v ■ V,.J((?'|A(xo)k")] {<i"\m\<i') dq"dq' . 

From we obtain 

2"(t;.V,J(g'|A(xo)k") = (g'|pSk">, 
where {q'\pQ\q") is given by ©. So far we have 

(^;- V,.J^(xo,t;^) = / {q'\pl\q"){q"\A{t)\q')dq"dq' = (I(t)),. (20) 



(21) 



Now, let A{x) be a vector function. Introduce 

A^(xo;;i) = lim \\n\\{v -V x„)A{xQ,t-h)\\ 

I— ►oo t 
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Using (|2U|1 we can also write 

A^{xo;h)=lim^ln\\{A{t)),.\\. (22) 
Taking limit ^ ^ in ^ we get 

A^(xo;0) = lim jln\\{vV,,)A{x{t,xom. (23) 

t— >oo t 

This formula is well defined for a wide class of vector functions A(x). To get 
the classical Lyapunov exponent one must make however the special choice 
A{x) = X. Then formula (|23jl becomes exactly 

In the quantum case symbol A{x) = x defines operator x which makes 
it possible to talk about the quantum trajectory ^ defined as the symbol of 
X{t) = U+it)xU{t) 

X{xo,t;h) = limTrX {t)pe = lim / X{x,t; h)We{x; xq)) dx . 

e^O e-+0 J 

Then formula H22() gives us definition of the quantum Lyapunov expo- 
nent, which coincides with ©, ((HI)- 

In view of (|19() lim^^^o Xixo,t; h) = x{t, xq), and the notion of quantum 
trajectory and (|2H) . (|22() allows also to write © in a form similar to 

A^ = A^{xo;h) = limlln{{vV^,)X{xo,t;h)}. (24) 

5. Radon transform and tomographic procedure. 

Let us make a few comments on the approach undertaken in [2] . To sim- 
plify exposition of some formulas and their geometrical meaning we assume 
in this section that n = 1, i.e q £ , x £ M? . 

Although the integral of the Wigner function for the density matrix 
is always 1, typically there are regions in phase space where W{x) is 
negative. This behavior results from the spectral expansion of the quantizer. 
Property 1 shows that the spectrum of A(x) is ±1 for all x. In detail 

A(x) = P+{x) - P_(x) , I = P+{x) + P_(x) (25) 

where P±{x) are the spectral projectors for the eigenvalues ±1. If x = 0, 
then A(0) is parity operator on Ti. and P±{0) are the corresponding even 
and odd projectors. When x ^ then P±{x) are the Heisenberg translates 
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of these operators, namely r(x/2)'^ P±(0)T(x/2). Applying to W{x) 
gives the Royer expansion 

wix) = l[\\p+{xM'-\\p.ixM'] . 

So whenever | |P+(x)'(/'| | < ||P_(a;)'(/'|| then W{x) is negative. This makes it 
impossible to interpret it as a classical statistical density. 

Nevertheless the two physically important projections of the W{x), namely 



/ 



Wiq,p) dp = , / Wiq,p) dq = 



are positive and form corresponding marginal distributions. Here ip{p) = 
J ^(g)g-w/^(f(^ is tJie wave function in the momentum representation. 
One may consider a family of projections with respect to all directions, 
in the phase space, not only the two given above. In this way we obtain 
Radon transform of W{x) 

RW{Q, tv) = j W{q, p) 6{Q -qi- pr]) dq dp . 

Here Q — q^ — p7] = is an equation of line in the phase space. Rewriting 
the above formula as 

RW{Q,tv) = — [ W{q,p)e~'''^^-'^^~P''Ukdqdp 
2vr J 

= ^ J dke~'^^ j dqdpW{q,p)e''^'"^^+''P'^^ 

one can see two remarkable things. First, the RHS can be viewed as a inverse 
Fourier of a characteristic function J W{q,p) exp{i{kq^ + kpr])}dqdp and 
thus is non-negative and represents a marginal distribution. Second, Radon 
transform is a composition of 2D-Fourier transform and ID-inverse Fourier 
transform, thus it is invertible. The inverse Radon transform can be written 
as 

W{q,p)= [ dQdidr]RW{Q,i,r])e'^^-'i^-P'^\ 



This means that knowing projections in all directions one can reconstruct 
Vl^(a;) from them by the inverse Radon transform. This property is widely 
used in computer tomography. 
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Due to invertibility of Radon transform the tomographic representation 
is equivalent to Weyl symbols representation. In particular the formula for 
the quantum mean takes form @ 



where A{^,r]) is inverse Fourier of A{q,p), and RW{Q,S,,7]) is the tomo- 
graphic distribution corresponding to W{x). 

Any of these representations finally lead to the same formula for the 
quantum Lyapunov exponent. 

This research was supported by grants to M.F.K. and T.A.O. from the 
Natural Sciences and Engineering Research Council of Canada. 
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